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Abstract 

We present a variant of the supereigenvalue model proposed before by Alvarez- 
Gaume, Itoyama, Manes, and Zadra. This model derives a set of three planar loop 
equations which takes the same form as the set of three anomalous Ward-Takahashi 
identities on the gaugino condensates recently derived by Cachazo, Douglas, Seiberg 
and Witten. Another model which implements J\f = 2 superVirasoro constraints is 
constructed for comparison. 
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Over the decades matrix models have contributed to our understanding of gauge fields and 
strings and their synthesis in different contexts. The recent proposal of Dijkgraaf and Vafa 
[1], (see [2] for prehistory to this), has enabled us to relate exact results on the superpotential 
and gauge couplings for a wide class of Af = 1 supersymmetric gauge theories with the planar 
limit of a matrix model in which the superpotential is taken as an ordinary potential. One 
may regard this relation as that connecting two faces of the notion of integrability of these 
theories in their infrared limit. The one face is seen as the integrability of Seiberg-Witten, 
i.e., Af = 2 theories [3] [4] with soft breaking superpotential as Whitham flows [5]. Here 
the aspects associated with gaugino condensates are not revealed immediately. The other 
face is, of course, the integrability of matrix models where the (planar) free energy provides 
an integral representation of the prepotential. Here, bosonic integrations of matrices have 
no way to capture the properties associated with supersymmetry and we have to input a 
given breaking pattern of the original gauge theory. It seems that we need to regard both 
gaugino condensates and vacuum values as moduli in order to gain unified understanding of 
this subject [6]. Studies of the Dijkgraaf- Vafa proposal from these perspectives have been 
continuing both on its structure and explicit calculations [8, 9, 6, 7]. For the other aspects on 
this subject, see an extensive list of references seen, for instance, in [6] and the last reference 
of [7]. 

The proposal of Dijkgraaf and Vafa has been put and consolidated in a field theoretic vein 
by Cachazo, Douglas, Seiberg and Witten [10] in the case of an adjoint chiral superfield and 
an arbitrary superpotential. In particular, they have derived an equation of the anomalous 
Ward-Takahashi identity associated with the Konishi anomaly [11] which takes the same 
form as the loop equation of the ordinary hermitian one-matrix model in the planar limit. 
In their discussion, there is actually a system of three Ward-Takahashi identities involving 
generating functions R(z),w a (z),T(z), which are all chiral operators related by constant 
fermionic shift of the gaugino field strengths. 1 The first operator R(z) realizes the proposal 
of [1] while the third operator T(z) is closely related to the Seiberg-Witten differential. 2 
The system of three Ward-Takahashi identities appears to offer a unified understanding of 
[1] and Af = 2 Seiberg-Witten theory. See [12] for a recent progress. In the ordinary bosonic 
one-matrix model, on the other hand, there is only one kind of resolvent whose Laplace 
transform is a loop operator, and is no hope within this model to obtain the second, which 
is fermionic, and the third resolvents carrying lower dimensions than the first one. A model 
which derives all of the three equations of [10] and which provides a constructive definition 
serving both for the planar processes and for the non-planar (and therefore) gravitational 
ones has remained elusive. 

1 This fermionic shift is allowed by the decoupling of the U(l) factor in the adjoint representation. 

2 This latter statement holds even in those cases in which unbroken nonabelian gauge symmetry survive. 
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To incorporate the counterpart of the fermionic shift operation, fermionic degrees of 
freedom are required. The supereigenvalue model proposed sometime ago in [13], which 
materializes discretized two-dimensional supergravity, (see [14] for some of the subsequent 
developments), has come to our attention. 3 In this letter, we will present a variant of 
the supereigenvalue model in which the ordinary grassmann Af = 2 superspace and its 
supercoordinates (A, 9, 9) 4 are introduced to label a set of eigenvalues \ and their grassmann 
partners 0*, 0j. We will show that this model captures the three equations of [10] as planar 
loop equations. 

Let us first consider the following partition function: 

z n= i ■■■ / II dXiddidSi I](A, - A, + 9 A - 9,9 jf exp (-- £ V(A l , i? 9,)) , (1) 

J J i=l i<j \ 9s i J 

where g s = -j| and 

V(Ai, 0i, 6i) = W^; g t ) + W (s \\ t ; Sift + 9 t W (s \X t] &) + W{\; g^Ofr (2) 

is an Af = 2 superfield in the terminology of the two-dimensional superconformal field theory 
in which the coefficients of the Taylor expansion in A« of the components W, W^ s \ W^ S \W 
are realized as the couplings ge, ^, gg. All of the singlet observables are obtained from this 
partition function by acting upon these couplings. We have fixed the normalization of the 
grassmann variable, so that the relative coefficient between \ — Xj and 9$i — 9j9j in the 
measure factor is 1. 5 

We bring this model into more direct contact to the proposal of [1] and its description as 
gauge theory in [10] by reducing this partition function. Let us set 

W (s) = W (s) = 0, W = W , (3) 

which is, of course, a nonsupersymmetric reduction. We regard the function W an n-th order 
polynomial in its argument. The model obtained this way is 

Z n= I • • • I II dXidOiffli H(\i - Xj + 9A - 9j9jf exp (-- ]T W(Xi + 9A)\ . (4) 

J J i=l i<j \ 9s i J 

3 Supermatrices have been known not to serve well in the old context of 2d supergravity. (See, for instance, 
[15] for this). For discussions in the present context, see [16]. There is a subtle issue of a decoupling of the 
superunitary group U(M\N) in the integrations over the supermatrix. These integrations in general do 
not reduce to those of the (super)eigenvalues. Furthermore, supertrace requires us to introduce oppositely 
charged particles in the Dyson gas description of the (super)eigenvalues. The supermatrices appear to us 
more appropriate to the quiver gauge theories and not to the pure Yang-Mills. 

4 6 and 9 are two independent grassmann variables and no complex conjugation is involved. 

5 In this letter, we set the value of the exponent of the measure factor to 2. For a general value /3, the 
coefficient of an order g s term in the loop equation below gets modified after a rescaling of the supcrpotcntial 
and the resolvent. 
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There is another supereigenvalue model which possesses Af = 2 supersymmetry, which 
we will present briefly in the end of this letter. In contract with the model (4), the model 
with Af = 2 supersymmetry is obtained by a different choice of the measure factor Yii<j(\ — 
Xj - (OiOj + 0i6j)) 2 . 6 

Let us briefly recall the logic which connects the Virasoro constraints at finite N with 
its loop equation in the ordinary hermitian one-matrix model [18]. Consider the partition 
function with all grassmann degrees of freedom removed in eq.(l). Setting up the system 
of Schwinger-Dyson equations associated with most general (polynomial) variations of the 
matrix is equivalent to inserting the Virasoro operators (which is a total derivative) iffl = 

N 

— X!dF^r +1 ) m > — 1 in the partition function and equating the resulting expressions to 

i=l 

zero. Summing over m by introducing an expansion parameter x, we obtain 

o = /-/nA(i:^)«p(^) • 

S%> EE -2£log(A ( -\-) + ilXA ( ) . 

i<j " s i 

Carrying out the differentiations and straightforward algebras, we obtain the finite N loop 

equation for the resolvent operator u^\x) = — fl^X/ir^A 7 : 

i=i 

JV j X — Ai 

Here ((■••)) denotes the averaging with respect to the partition function. The absence of the 
g s ((ujft {%))) term is related to c = 1 nature of the Virasoro constraints. In the planar limit 
in which the correlators of the singlet operators factorize, we obtain l)( b \x) = lim u^\x): 

N->oo 

J B \x) 2 + W'(x)^ B \x)+ \im((iy j W ' {x) ~ W ' {Xi) }} = . 

We now apply the above procedure to our model given by eq. (4). We find that the 
appropriate realization of the Virasoro operators which involve the grassmann eigenvalues 

6 These two alternative choices may be characterized as follows: by requiring the antisymmetry in 
and Af = 2 neutralness, possible bilinear combinations are either OiOi — 9j9j or 9i0j + 9{9j. The difference 
is that the former is odd under the exchange of 9 and 9 while the latter is even. Thus the former may be 
called chiral choice and the latter non-chiral. In fact the combinations \ ±99 are the natural coordinates for 
chiral and anti-chiral superfields respectively. 
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as well is 

c - -E^r- +7(m+ i)i:(-| Si+ -lft)Ar, ™>-i 

«*) - - E + 7 E ( A 9j + ^) ^ . ,5) 

This form can be motivated by realizing the M = 2 superconformal algebra in terms of 
(super-)differential operators. The coefficient 7 is then determined: 

7=5. (6) 

Inserting — which is a total derivative, in Z^, carrying out algebras, and equating 
the resulting expression to zero, we obtain 

=d-»r)«E(^» 

- l«£ (jh^A: + "J^r + + ¥i))) (7) 

+ 2«E + 7^ + glog (A, - A, + a - fcfc 

The first line of this equation vanishes with our choice 7 = |. The second line with this 
choice is 

1 _ if' (A, + OjOj) 

The calculation of the third line is a little less trivial to see. After antisymmetrization and 
isolation of purely bosonic term, we obtain 

//y* \ u _ //y ^i^i ~ = \\ 

U ^ (x - A,)(x - A,) /; U ^ (x - A,)(x - A*) (A,- - A, + - QAY 

1 o-vZ/V f ^? e Jh_ \ 1 _ u 

m £J \{x - A,) 2 (x - \Y) \j - A, + 0^ - QSi ' 

The cancellation of the 4j contribution from the second and the third terms of this equation 
in the asymptotic expansion again requires 7 = \. With this choice, the third line of eq.(7) 
is 

% 3 ( X - Xj - e&Xx - \i - ea^ ■ (8) 



-5- 



N 



Introducing the resolvent operator u>ft(x) = ~9s^2 ^z\^Zgj. i we obtain 



= ((^(x) 2 )) - g.«a^(x)» + + ((| £ ^ f (A ; t 

= • (9) 

The fourth term is a polynomial in x. In the planar limit u>(x) = limcj^rr), this equation 
reduces to 

^(x) 2 + W(x)o;(x) + hm((-^ ^ 1 ^ x K nn ' )) = • (10) 

iv j a; — Aj — 

So far we have managed to obtain a single loop equation for the resolvent. The observables 
of our model are, on the other hand, 

i i i i 

With the single expansion parameter x, we are not able to generate all of these which are 
independent. In order to overcome this point, we insert 

following the insertion of the Virasoro operators in eq.(7) and carry out the remaining pro- 
cedure. The grassmann expansion parameters 77, fj have now been introduced. Equivalently, 
we can make a fermionic shift of the grassmann eigenvalues in eq.(9): 

Both the resolvent operator oj^{x) and 

s cS = - £ w(\ t + eA) - 2 J2 io g (A, - \j + e t e t - ofa) (11) 

9s i i<j 

undergo this shift while d(9i + rj) — d$i , d($i + fj) = dO^. In either way, we obtain another 
set of Schwinger-Dyson equations as equations multiplying 77, 77, 777/; 

= ((C 2 (x)+Sr h C 1 (x))) , 

= ((AW+5f b A(x))) , (12) 

= ((jc 3 (x) + (sit+sr b sr h -((si» h sr h )^jc 1 (x)+ ^s^+s^^x))) , 
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where 



C 2 (x) = 2ipx(x)u>x(x) - g,il>'jf(x) + W'(x)ip^(x) + p^(x) , 

Czix) = 2'4) f j{x)ujf { {x) - g s i>'fj{x) + W'{x)i)^{x) + p^(x) , (13) 
£ 3 (x) = 2t^(x)u^(x) + 2i)fj(x)i)^{x) - g s t'fr(x) + W'{x)t^(x) + c^(x) , 



and 



N * (x-\-BA) N * (x-Xi-OiOi) 

S W^(x; Xj + 6 A) g ^ - Xj + gA) 

N * (x-x-oa) n * (x-Xi-eA) 

We have introduced 

W^ +1 \x- Aj + 6^) ee - £ (X " A< " W ^ (m) ^ + W • ( 15 ) 

We see that the quantities seen in eq.(14) are polynomials in x. Another bosonic resolvent 
and two fermionic resolvents are defined respectively by 

/ 

t#(x) ee -g s J2 



I 20^ 



(x-x t -eA) 2 (x-K-eA) 3 ' 



^(x) EE -gsJ2~, ^T2' M x ) = ~9a^2 



(x-x t -eA) -^{x-K-eAf 



In the planar limit, the terms involving the singlet operators coming from the variation 
of the action which we have denoted by 

*-s:^. da) 

factorize and become ignorable. We obtain 

2^(x)u(x) + W\x)ip(x) + \im((px(x))) = , 

2ip(x)u(x) + W\x)i>(x) + lim((p#(a;))) = , (17) 
2t(x)u(x) + 24>(x)ip(x) + W'(x)t(x) + lim((c J v(x))) = . 

As is announced, the system of three planar loop obtained in eq.(17) takes the same form as 
that of the equations of Konishi anomaly in [10]. 



So far, we have dealt with the case in which the underlying gauge group is U (N); this has 
led us to the supereigenvalue generalization of the hermitian one-matrix model. It is easy to 
generalize to the cases in which the underlying gauge group is SO(2N)/ Sp(2N), SO(2N+l). 
(See the recent discussion [17]). The measure factor is changed into 

N 

Y[(\i - Xj + 9 l 9 l - 9 1 9 1 ) 2 (X t + Xj + OA + Q 3 Qjf , for SO(2N) 

i<j 

N N 

JJ( A . - Xj + OA - 9,9 ,) 2 (A, + Xj + 9 t 9 t + 9j9j) 2 U( X i + e M > for M™), SO(2N + 1) 

i<j i 

These cases can be handled by changing the set which integer labels i,j have belonged to 
into the new ones: 

i,je{l,2,---,N} => i,je{-7V,...,-2,-l,l,2,...,JV}, iorSO(2N) 
i,je{l,2,-..,N} i,je{-JV,...,-2,-l,0,l,2,...,iV}, iorSp(2N),SO(2N) 

and by demanding 

X-i = —Xi , Ao = . 

The summations and the products in the previous equations are now taken with respect 
the integers belonging to these new sets. The potential effectively becomes a polynomial 
consisting only of even powers. The remainder of the discussions are the same. 

Finally, let us discuss the supereigenvalue model which implements Af = 2 super Virasoro 
constraints. The action is given by the M = 2 superfield (2). In the partition function (1), 
the measure factor is replaced by 

H(X t - X, - (9 t 9, + 9 t 9j)) 2 , (18) 

i<j 

where each factor is the super-distance invariant under the supertranslation generated by the 
supercovariant derivatives. The (super)-differential operators to be inserted in the partition 
function are 
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and £(x) (the case 7 = |) defined in eq. (5). They can be organized into a current superfield. 
The same is true for the resolvents of this model: 

n(l ,rf,- S ,E ^M . (20) 

j X A j 

Introducing the supercovariant derivatives D = ^- + r}-^ , D = J= + 77^, we have been 
able to derive a manifestly M = 2 supersymmetric loop equation. Space only permits us to 
present its final form: 

= ((DQDQ)) - g s ((Q')) + (([dVDQ + DVDQ))) + ((F)) . (21) 

Here F(x, 77,77) is a polynomial in x given explicitly in terms of the components of V. Fac- 
torization of the two-point functions in the large N limit (g s — > 0) gives us an M = 2 planar 
loop equation. 

We plan to give full details of the contents of this letter together with other issues in a 
future publication. We thank Alyosha Morozov for helpful discussion on this subject. 
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